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The purple bronze Lio.gMo60i7 is of interest due to its quasi-one-dimensional electronic structure 
and the possible Luttinger liquid behavior resulting from it. For sufficiently low temperatures, it is 
a superconductor with a pairing symmetry that is still to be determined. To shed light on this issue, 
we analyze a minimal Hubbard model for this material involving four molybdenum orbitals per unit 
cell near quarter filling, using asymptotically exact perturbative renormalization group methods. 
We find that spin triplet odd-parity superconductivity is the dominant instability. Approximate 
nesting properties of the two quasi-one-dimensional Fermi surfaces enhance certain second-order 
processes, which play crucial roles in determining the structure of the pairing gap. Notably, we find 
that the gap has more sign changes than required by the point-group symmetry. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27.+a, 74.20.Rp 


I. INTRODUCTION 

Understanding unconventional superconductors, char¬ 
acterized by a sign-changing order parameter, remains 
a challenge 0 , 1 ]. Most known unconventional super¬ 
conductors (cuprates, organic charge transfer salts, iron 
pnictides) involve pairing in the spin-singlet channel. The 
leading candidates for spin-triplet pairing are strontium 
ruthenate [3] and the heavy-fermion compound UPt 3 - In 
general, it remains unclear how the pairing symmetry is 
determined by the interplay of strong electron correla¬ 
tions, low dimensionality, spin fluctuations 0 , multiple 
bands Q, Hund’s rule coupling, band structure effects 
[fjj, and spin correlations in proximate Mott-insulating 
states 0 . 

Lio.gMoe0 17 , known more commonly as a “purple 
bronze,” is a layered transition metal oxide that ex¬ 
hibits several exotic properties arising from the effec¬ 
tive low dimensionality and strong electron correlations. 
At high temperatures, the dynamics are primarily one¬ 
dimensional and the system exhibits some properties that 
have been interpreted in a Luttinger liquid framework 0 
0 . At lower temperatures the metallic phase is difficult 
to characterize, possibly due to incipient charge density 
wave order [ll],[l2|. Still, for magnetic fields perpendicu¬ 
lar to the easy transport axis (the crystalline b axis), the 
magnetoresistance grows without bound as a function of 
the field strength, as consistent with an open Fermi sur¬ 
face 0 - There are several unconventional aspects of 
the superconducting state that are worth noting. First, 
the transition temperature T c decreases with increasing 
disorder (residual resistivity) 0 , a signature of non-s- 
wave pairing (for a review, see Ref. QjJ). Second, the 
upper critical field exhibits a great deal of anisotropy, 
consistent with a quasi-one-dimensional system. In par¬ 
ticular, for a magnetic field applied along the b axis, the 
upper critical field exceeds the Chandrasekhar-Clogston 


paramagnetic limit 0, suggesting that the supercon¬ 
ductivity may be spin-triplet in character. However, this 
requires ruling out alternative explanations such as the 
effects of spin-orbit coupling and quasi-one-dimensional 
fluctuations 0 - 

A well-studied example of unconventional pairing in 
quasi-one-dimensional systems is the organic Bechgaard 
salts (TMTSF) 2 X. In these materials, there is evidence 
in favor of nodal, spin-singlet pairing [l 8 j . However, in 
the presence of a strong in-plane magnetic field applied 
along the chains, superconductivity in these systems also 
appears to exceed the Chandrasekhar-Clogston limit. If 
a modest magnetic held can induce a transition of this 
sort, it is quite natural to infer that there is a near de¬ 
generacy of pairing in the spin-singlet and spin-triplet 
channels. Indeed, weak-coupling considerations based on 
Fermi surface nesting do indicate such a near degeneracy 
in these systems [lj| . Whether the pairing is singlet or 
triplet remains a subject of debate [201 ]. Here we explore 
the role of Fermi surface nesting, multiband effects, and 
low dimensionality for unconventional superconductivity 
in the context of LionMoeOir. 

Despite evidence of strong correlation effects in purple 
bronze, we take here a weak-coupling perspective. The 
advantage of such an approach is that the solutions ob¬ 
tained are well controlled, and robust trends suggesting a 
certain pairing symmetry can easily be identified in this 
limit. The qualitative information obtained this way may 
be relevant to the experimentally relevant intermediate¬ 
coupling regime. Indeed, for the cuprates it has been 
found that the weak-coupling predictions of d- wave su¬ 
perconductivity are consistent with spin-fluctuation ap¬ 
proaches and strong-coupling approaches such as RVB 
theory [d]. For simplicity, we consider a four-band model 
Hamiltonian with on-site repulsive electron interactions. 

Our principal findings may be summarized as follows: 
Unlike the case of Bechgaard salts, the near degeneracy 
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FIG. 1. Tight-binding lattice model. Each circle corresponds 
to a single Mo atom and there are four atoms per unit cell. 
Filled and empty circles denote two types of crystallograph- 
ically inequivalent Mo atoms, i.e., Mo(l) and Mo(4), within 
a layer of Lio.gMo60i7. The horizontal orange line defines a 
plane of mirror symmetry, and the green dot an axis of C 2 - 
rotation symmetry. (The reflection is only an approximate 
symmetry in the real material, which is slightly monoclinic 

m) 

between the dominant triplet and the dominant singlet 
solutions is absent; the trivial pairing channel is signifi¬ 


cantly penalized, as the bare on-site repulsion takes effect 
in multiple sublattices. Significantly, the leading pair¬ 
ing instability is in the odd-parity channel; there is a 
sign change in the gap function (not required by sym¬ 
metry) across the two nearly degenerate Fermi surfaces, 
and additional accidental nodes develop in the inner sur¬ 
face when the parameter inducing the splitting of the 
two Fermi surfaces as well as their warping is sufficiently 
large. 


II. THE MODEL 

In Li 0 .gMo 6 Oi 7 , the low-energy electronic degrees of 
freedom reside on ladders built from d xy orbitals of Mo 
atoms [2^ - [24| . which run along the crystallographic b 
axis. Based on this observation, the electronic struc¬ 
ture was recently described by a four-band tight-binding 
model near quarter filling, i.e., n e i = 1.9 out of 8 per 
unit cell, in Ref. EH (Similar, but slightly different mod¬ 
els have been presented in Refs. EqI and [25].) In this 
model, bonds within each chain are formed by the hop¬ 
ping amplitude t , whereas t 1 - connects two nearest chains 
to make a ladder, and neighboring ladders are coupled by 
t' (Fig. [TJ. The resulting single-particle Hamiltonian can 
be represented in the momentum space as follows: 


hr A k) 


( 


0 


t'e ikc (1 + e ikb ) 
\ t( l + e ikb ) 


0 

f(l + e ikb ) 

0 


t'e~ ikc {l + e~ ikb ) t{ 1 + e~ ikb ) \ 
t(l + e~ ikb ) 0 

0 t 1 - 

t 1 - 0 / 


(1) 


Here, r and t' are sublattice indices, and in the basis used 
to construct the above matrix, the order of sublattices 
is 1, 4, 1', and 4'. (See Fig. [T]) We denote the n th 
eigenvalue of /i rT /(k) as a T (n, k) and the corresponding 
band dispersion as e(n, k). 

So long as |i J -|,|t'| <C |t| and the filling fraction is 
near one-quarter, two nearly degenerate Fermi surfaces 
exist. They are open along the c axis, and there is an 
approximate interband nesting betweeen them at Q = 
(0,7m e i/2) = (0, 0.957t) (Fig. [2]). To identify a robust 
tendency that does not involve a fine-tuning of the tight- 
binding band structure, we introduce the warping param¬ 
eter rj w , where t ± = —0.04:8r] w t and t' = 0.072 r] w t. As rj w 
is increased, the splitting between the two Fermi surfaces 
and the warping of them become larger; at the same time, 
the nesting becomes less perfect. The band structure is 
reduced to that described in Ref. Ell at n,„ = 1. 

We further add on-site repulsive interactions that lead 
to a Hubbard model with the following Hamiltonian: 

H = H 0 + Hint, 


^ ^ (k) Cj crcr Ci cr / fT , (2) 

k ,t,t' ,<j 


H lilt — jy ^2 c kiTt C k2Tl Ck 4-G c k3-rT) (3) 

with the understanding that k 4 = ki + k 2 — kg. Spin 
indices are denoted as cr and cr', and number of unit cells 
as N. 

The weak-coupling regime corresponds to the case that 
|t|, 1^1, |t'| U. In this limit, it is desirable to work in 
the band basis that makes Hq diagonal as follows: 

C]ct<t — \ a T {n, k) c n k(j? (4) 

n 

and 

H o = e ( n ’ k ) c ika c «k ff - 

n,k,tr 


( 5 ) 
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III. PERTURBATIVE RENORMALIZATION la Kohn and Luttinger [27} . The effective interaction at 

GROUP (RG) METHOD energy scales asymptotically close to the Fermi surface is 

given by the following expression: 

We implement a renormalization group (RG) method 
Haim to investigate superconducting instabilities a 

I 


r n (n, k;n',k ) = r^(n, k;n',k ) + rj^(n, k; n', k ), 
r^ ) (n,k;n , ,k ) =U^2 |a r (n,k)| * 2 |a r (n',k )| 2 , 




k; n', k ) = 



/(e(m', p) - n) - } (e(m. p + k + k ) - /x) 
e(m', p) - e(m, p + k + k ) 


-(n, k)* a T [m , p + k + k ) a T (m', p)* a T (n r , — k ) 


( 6 ) 


Here, / is the Fermi function at zero temperature, and n 
is the chemical potential. It represents a process where 
two electrons of opposite spins scatter from (n', ±k ) to 
(n, ±k) without spin flip. It is unnecessary to consider 
the corresponding object for equal-spin electrons: due to 
spin-rotation invariance, Eq. ([6]) accounts for both the 
singlet and triplet channels. 

Different pairing channels are identified by solving the 
following integral eigenequation defined along the Fermi 
surface: 


E 


- ds( ' kn ’\ f T n (n , k n ; ri ,k„,) i/} a (ri, k ra ,) 

(27r) 2 v(n', k„,) (7) 

— AcK^a k n ) . 


Here, k n denotes a point on the Fermi surface in the n th 
band, ds( k„) is an infinitesimal length element of this 
Fermi surface, and v is the Fermi velocity. The pair wave 
function tjj a is classified according to its parities under 
the transformations k c —> —k c and kb —> —kb■ These 
symmetry operations in the momentum space are derived 
from the symmetry in the real space shown in Fig. [T] 

Negative eigenvalues of indicate the sectors in 
which the effective interaction is attractive. These at¬ 
tractive interactions grow under further renormalization. 
The eigenfunction with the most negative eigenvalue (de¬ 
noted as Ao) is responsible for superconductivity; pairing 
instability sets in at the energy scale of We -1/1 1 A ° I, where 
W = 4t is the bandwidth. This energy scale is also iden¬ 
tified with the transition temperature T c , and the gap 
structure is inherited from i/jq (n,k„). 

Notice that in the asymptotic weak-coupling limit, the 
bare on-site repulsion in Eq. © is infinitely larger 
( 2 ) 

than the other term . Hence, attraction can arise only 
if the pair wave function completely avoids the effect of 


-.(i) 

- ft ’ 


i.e., 


belongs to its null space: 


E 


ds( k„) 


(2n) 2 v(n, k„) 


\a T (n, k„)|V a (n,k„) = 0 (8) 


for all four possible values of the sublattice index r. In 
practice, one should solve Eq. © with replaced by 

( 2 ) 

ri/ while enforcing the above conditions. 

It should be noted that Eq. (j8]) is automatically satis¬ 
fied by pair wave functions that are odd under k c —»• —k c 
or kb —> —kb- To see this, observe that |a T (n,k„)| 2 is an 



k c /n 

FIG. 2. The Fermi surface and the sign structure of the 
(triplet) gap function for rj w — 1. The inner and outer 
curves correspond to Fermi surface sheets associated with 
distinct bands. The vertical arrow depicts the nesting vec¬ 
tor Q = (0, 0.957r). The oblique arrows denote an example of 
two points on the Fermi surface whose crystal momenta ap¬ 
proximately sum to Q; the gap function changes sign between 
these points although not required by symmetry. 
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even function of both k c and kb. This fact is a result of 
the reflection (see Fig. © and time-reversal symmetries. 
When ip a {n , k„) is an odd function of either k c or kb , the 
contributions from different quadrants of the Brillouin 
zone cancel out in the integral in Eq. (|8j>. On the other 
hand, the four conditions given by Eq. © lead to addi¬ 
tional constraints on solutions that are even in both k c 
and kb- However, due to the ^-rotation symmetry, Eq. 
© becomes identical for r = 1,1' as well as for r = 4,4' 
(Fig. [T]), and the number of constraints is two as opposed 
to four. 

In addition to constraints set by the bare repulsion, 
nesting properties of the Fermi surface play an impor¬ 
tant role in determining the gap structure. In our case, 

Ti / (n, k; n', k ) in Eq. © is considerably enhanced when 

- 

the momentum transfer k+k is close to the approximate 
nesting vector Q = (0.957 t, 0) (Fig. ©). Moreover, from 

( 2 ) 

its definition, we see that TW is always positive. Sum- 

ming up, the nesting condition k+k ss Q leads to a large 

- / 

positive off-diagonal matrix element connecting k and k . 
Therefore, one can expect that a pair wave function with 
a large negative eigenvalue tends to change sign between 
two points on the Fermi surface whose crystal momenta 
sum to Q. 

Together with the repulsive Hubbard interaction, the 
nesting property of the Fermi surface shown in Fig. [2] 
causes spin density wave (SDW) fluctuations of wave vec¬ 
tor Q. In view of this fact, the previous paragraph is 
about how SDW fluctuations enhance superconductivity 
and determine the gap structure. Notice that within our 
weak-coupling approach, such nonsuperconducting fluc¬ 
tuations are never strong enough to compete with su¬ 
perconductivity unless an extreme fine-tuning is made. 
That is, as long as the band dispersion is invariant un¬ 
der k -+ —k, the perfect nesting in the particle-particle 
channel makes superconductivity dominant over other 
types of instabilities except at certain special points in 
the parameter space (e.g., Van Hove singularity, per¬ 
fect nesting in the particle-hole channel). The one¬ 
dimensional limit (r\ w = 0) corresponds to such a special 
point, but we avoided it. Then, as we know that the 
ground state is superconducting, it is indeed a valid ap¬ 
proach to investigate the form of pairing that leads to the 
strongest instability of the Fermi sea without considering 
nonsuperconducting instabilities. Correct description of 
the ID physics at rj w = 0 is beyond the scope of this 
study. 


IV. RESULTS OF THE WEAK-COUPLING 
ANALYSIS 

Figure [3] shows the dominant pairing strength in the 
triplet and singlet channels as functions of the warping 
parameter r ] w . (We have plotted the dimensionless quan¬ 
tity A = A jjT-) To obtain this result, we numerically 
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FIG. 3. Dominance of triplet pairing is insensitive to the 
extent of warping of the Fermi surface. The pairing strength 
in the triplet and singlet channels is shown as a function of the 
parameter ri w that determines warping of the Fermi surface. 
Reference Hll estimates that r/ w ~ 1. 


solved Eq. © by discretizing each of the two Fermi 

( 2 ) 

surfaces into 128 segments and evaluated Ti/ in Eq. 
© by discretizing the Brillouin zone into a 2048x2048 
grid. f28l] The primary conclusion to draw from the weak- 
coupling analysis is that the pairing instability is in the 
spin-triplet channel. For a broad range of r) w , the domi¬ 
nant pairing solution is a “p^-wave state” that is even in 
k c and odd in kb . The associated pair wave function for 
Vw = 1-0 is shown in Fig. [I] 

Notice that there are more sign changes in the pairing 
gap than dictated by parities under k c —>• — k c or kb —> 
—kb- Most importantly, it has opposite signs between 
two nearly touching Fermi surface sheets. This fact is 
consistent with the observation in the previous section 
that the pair wave function has a propensity to change 

.w ~ 

sign between two points k and k satisfying k + k « Q 
(Fig. ©. 

Another notable feature is the gap suppression near 
k c = 0, as seen in Fig. [I] A gap minimum exists on the 
outer Fermi surface at k c = 0 for the entire range of r/ w 
shown in Fig. [3] and also on the inner Fermi surface for 
Vw 0-8- For larger warping, two closely spaced acci¬ 
dental nodes develop around k c = 0 on the inner Fermi 
surface. The reason for this gap suppression can be again 
understood from the nesting condition k + k ss Q. From 
Fig. m one sees that most of crystal momentum pairs sat¬ 
isfying the nesting condition exist across the inner and 
outer Fermi surfaces. However, near k c = 0, the two 
Fermi surfaces almost coincide, and a pair of crystal mo¬ 
menta on the same Fermi surface, as well as a pair across 
different Fermi surfaces, can meet the nesting condition. 
Hence, there is frustration in the resulting sign change of 
the gap, leading to the suppression. 

We remark that from the above consideration of the 


Triplet 

Singlet 
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FIG. 4. The pair wave function for rj w = 1 on the inner and 
outer Fermi surfaces in the upper half of the Brillouin zone. 


nesting properties, one may expect that that there exists 
a closely competing singlet solution, which has the same 
gap structure in the upper (or lower) half of the Brillouin 
zone as the dominant triplet solution but is even in both 
k c and kb- However, the numerical result shows that this 
is not the case. Indeed, for solutions that are even in 
both k c and kb , one has to additionally apply constraints 
given by Eq. ©. As a result, what appears to be a 
nearly degenerate singlet solution is in fact significantly 
disfavored. 


V. ROLE OF LONGER-RANGE COULOMB 
INTERACTIONS AND SPIN-ORBIT COUPLING 

We also briefly considered the possible effects of off¬ 
site Coulomb interaction parameters introduced in Ref. 
fill , and set them to be of order U 2 /t [Hj]. However, if 
one neglects screening, there are a large number of these 
parameters. A representative investigation of the full 
parameter space would require another separate study. 
Here, we just note that we did find that these interactions 
could change the pairing symmetry and/or the number 
of accidental nodes in the pairing amplitude. We also 
found that the repulsion between diagonal neighbors has 
a much more significant effect than that between nearest 
neighbors along or in the direction perpendicular to the 
zigzag chain. Whether this sensitivity to longerOrange 
interactions is an artifact of the weak-coupling analysis 
may only be revealed by numerical investigations of the 
intermediate- to strong-coupling regime. 

Another important question is the direction of the d 


vector associated with the triplet pairing. This will be 
determined by spin-orbit coupling and we plan to address 
this in a future study. 

VI. PROPOSED FUTURE EXPERIMENTS 

One signature of triplet pairing is the temperature in¬ 
dependence of the Knight shift in nuclear magnetic res¬ 
onance (NMR) below T c [3:]. Hence, we propose mea¬ 
surements using 7 Li, 17 0, and/or 95 Mo NMR. Recently, 
7 Li NMR measurements were performed on the normal 
state of this material i[30]. In the past, 95 Mo NMR has 
been difficult due to the very low resonance frequency. 
However, recent advances in high-field NMR have seen 
its use in studies of a range of molybdenum oxides [3lj ]. 
Evidence for the nodes in the superconducting energy 
gap will be a power-law temperature dependence of ther¬ 
modynamic quantities such as the specific heat and the 
superfluid density at temperatures much less than T c . 
A recent theoretical study investigated the difference in 
the temperature dependence of the upper critical field for 
triplet states in a quasi-one-dimensional superconductor 
with and without nodes j32|. However, the differences are 
small and subtle, making definitive conclusions about the 
presence or absence of nodes difficult. 

VII. CONCLUSIONS 

We used asymptotically exact perturbative renormal¬ 
ization group methods to analyze a minimal Hubbard 
model for Lio.gMogOiy. We found that spin-triplet odd- 
parity superconductivity is the dominant instability, with 
a gap having more sign changes than required by the 
point-group symmetry. Hopefully, our study will stim¬ 
ulate further experimental studies of the superconduct¬ 
ing state of this material, to find definitive evidence for 
triplet pairing and/or nodes in the energy gap. 

Note added. Recently we learned of a work by Lera and 
Alvarez j33{ that performed a multi-orbital RPA analysis 
of a similar Hubbard model and reached similar conclu¬ 
sions to ours. 
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